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Abstract 

We present a variational study of employing the trigonometric basis functions satisfying peri- 
odic boundary condition for the accurate calculation of eigenvalues and eigenfunctions of quartic 
double-well oscillators. Contrary to usual Dirichlet boundary condition, imposing periodic bound- 
ary condition on the basis functions results in the existence of an inflection point with vanishing 
curvature in the graph of the energy versus the domain of the variable. We show that this boundary 
condition results in a higher accuracy in comparison to Dirichlet boundary condition. This is due 
to the fact that the periodic basis functions are not necessarily forced to vanish at the boundaries 
and can properly fit themselves to the exact solutions. 
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1 Introduction 



The phenomena of quantum tunneling in the double-well potential (Fig. [T]) 

V(x) = -kx 2 + Xx 4 (fc,A>0) (1) 

and finding the energy eigenvalues and its stationary states is a long-standing and well-known problem in 
quantum mechanics. The interest in this problem ranges from various branches of physics to chemistry. 
After the seminal works by Bender and Wu [T], numerous papers have been appeared in the literature 
to tackle this problem. In spite of its simple form, this problem presents a non-trivial model in both 
quantum mechanics and quantum field theory. In fact, this model can be considered as 0+ 1-dimensional 
quantum field theory X(j) which has no spatial dimensions. 
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In chemistry, this type of potential corresponds to two equilibrium positions or molecular config- 
urations. The periodic inversion of ammonia (NH3) is a well-known example of quantum tunneling. 
The internal rotation in CH3CH3 molecule from one configuration to another is a good example of a 
symmetric double- well potential. On the other hand, the asymmetric H2O2 hindered rotor or beryl- 
lium dicyclopentadienyl represent asymmetric double well potentials. Ammonia (NH3), cyanamide 
(NH2CN), PH3, and AsHs are chemical examples of symmetrical double well potentials where atoms 
can tunnel through the barrier. 

It is obvious that it is not possible to set up a manageable perturbation theory even for an infinitesimal 
value of the coupling strength A to calculate the exponentially small energy splitting between quasi- 
degenerate levels in a double-well. This stems from its asymptotic divergent nature [2] due to the 
eventual dominance of the perturbative correction over the unperturbed contribution for large amplitude 
of oscillation [T] . This has been the driving force behind the development of non-perturbative methods for 
this kind of problems. Among them are the semiclassical approximation [3J, finite-difference technique [3], 
hypervirial recurrencerelation scheme [5], renormalized hypervirial-Pade scheme [6], variational matrix 
solution [7], Rayleigh-Ritz variational method supplemented by the Lanczos algorithm 8 , Brillouin- 
Wigner perturbation theory based on shifted oscillator variational functions [9], instanton method [TOl 
111] , transfer matrix method |12j , uniform approximation of the logarithmic derivative of the ground state 
eigenfunction [13] . and many other specific methods [TH [T51 [TH1 [T7] . For the double- well potential, even 
quite sophisticated techniques, such as the renormalized hypervirial-Pade scheme [6], fail to calculate 
the semi-degenerate low energy states. 

The WKB approximation is widely used for its simple mathematical form, but the results are known 
to be inaccurate due to its inherent defect in the connection formula. One can take the quadratic 
connection formula instead of those related to the Airy functions to modify the WKB result for the 
ground state [TTJ [IB] . Some other refinements have been developed to improve the accuracy of WKB 
by changing the phase loss at the classical turning points [HI [20] . The anharmonicity is also taken into 
consideration in the case of small separation distance between the wells [20 . 
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In the instanton method, one uses imaginary time path integral to obtain the classical action the so- 
called Euclidean action. Qualitatively, this method is useful for understanding the quantum tunneling, 
which has no classical counterpart. However, the quantitative calculation of the energy splitting in the 
double-well potential by this method is inaccurate because the Euclidean propagator can be obtained 
only in the limit of infinite separation between the two potential minima |10j , which corresponds to zero 
tunneling probability. Thus, the validity of the instanton approach is restricted to the case in which the 
two potential minima are far apart, and its accuracy is expected to be reduced as they become close. 

Various variational methods are usually used in quantum mechanics, statistical mechanics and field 
theory In quantum mechanics, variational parameters are incorporated into trial wave functions and trial 
Hamiltonians. The Rayleigh-Ritz method is the minimization of the ground state energy with respect 
to these variational parameters. The applications of this formalism using various set of basis such as 
harmonic-oscillator [3T] , Chebyshev polynomials [35] , hypervirial theorems [53] , and the coherent-states 
[24) . have been already appeared in the literature. 

In Refs. [25j [26l [27l [28] , it was shown by numerical results that the eigenfunctions which obey Dirichlet 
boundary condition 

*(-L) = *(L) = 0, (2) 

can be effectively used to find the spectrum of an unbounded problem. Then, a critical distance L was 
defined, and it was shown that the low lying energy levels E n (L) are equal to those of L — oo with high 
accuracy, if the boundedness parameter L is in near vicinity of L. In these works, they have employed 
the Rayleigh-Ritz variational method and have chosen the trigonometric functions (particle-in-a-box) 
as basis set, in order to solve double- well potential and some other anharmonic symmetric oscillators 
25, 26, 27, 28 . They have shown that the trigonometric functions are a suitable basis-set in variational 
calculations in order to obtain highly accurate results. Moreover, the model is simple, fast-convergent, 
efficient and works for different kind of potentials [25, 26, 27, 28 . 
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In this paper, we will also use trigonometric basis functions, but with periodic boundary condition 



In this case, contrary to Dirichlet boundary condition ([2]), E n (L) has not any relevant minimum value. 
In fact, for periodic boundary condition E n (Q) is zero, whereas E n (0) is infinite for Dirichlet boundary 
condition. Note that, imposing periodic boundary condition on e.g. the ground state wave function, 
makes it more accurate than imposing Dirichlet boundary condition. This is due to the fact that the 
Dirichlet boundary condition enforces the wave function to vanish at x = ±L, but the periodic boundary 
condition lets the wave function to fit itself more closer to the exact solution which is not necessarily 
zero at the boundaries. Therefore, we show that using periodic boundary condition, the same accuracy 
can be obtained with the smaller number of basis. 

The paper is organized as follows: In Sec. [2 we outline the general calculation scheme using the 
periodic boundary condition. In Sec. 02 we apply this method to the case of double-well potential and 
compare our results with the previous near-exact ones. For the case k = 1, we find the accurate energy 
eigenvalues for various value of A and obtain the optimum value of L in each case. Also, we outline 
the scaling properties of the Hamiltonian to generalize the results to more complicated forms of the 
potential. In Sec. Efl we state our conclusions. 

2 The variational method 

Let us consider the following dimensionless time-independent one-dimensional Schrodinger equation 



where V(x) = V(—x). Since the potential is an even function of x, to avoid large matrices, we can use 
even and odd basis functions separately 



= #(£). 



(3) 



' d 2 i 



(4) 




m = 0, 1, ... 



(5) 
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Figure 1: One-dimensional double-well potential where k = 1 and A = 0.03. The ± signs on each energy 
level refers to the parity of the corresponding eigenstate. Starting from the ground state up, the energy 
splitting between the levels below the cusp are approximately {4.3 x 10~ 6 , 6.1 x 10~ 4 , 3.0 x 10~ 2 }. 



and 
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m = 1,2, 



(6) 



Now, we can expand the solution using these orthonormal basis sets 



= ^^(x), (7) 
m 

where i — 0, 1 denotes even and odd solutions, respectively. Also, we have the following expansion 



(8) 



where B m are coefficients that can be determined once V(x) is specified. By substituting Eqs. (0 and 
© into Eq. (0} and using the linear independence of the basis functions, we have 



frmr\- 



Am, ~\- Bm. — E A r . 



(9) 



Moreover, we can find B m 's using Eq. (|SJ) and (J7J, namely 



where 



C m ,m' = I 4>m{x) V(x) 4> m '(x) dx. 



(10) 



(11) 
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Therefore we can rewrite Eq. © as 




A-m ~l~ ^ ^ Cmj-fi' A in ' — E A. 



m ■ 



(12) 



By selecting a finite subset of the basis functions, i.e. choosing the first N, Eq. (fT2j) can be written as 



Therefore, the eigenvalues and eigenfunctions of the Schrodinger equation (j4|) are approximately equal 
to the corresponding quantities of the matrix D. 

3 Results 

To establish the method, first let us study an exactly solvable case where k = — 1 and A = 0. Simple 
Harmonic Oscillator has a well-known energy spectrum E n = 2n + 1. Figure [2] shows its ground state 
energy versus L for both periodic and Dirichlct boundary conditions. As apparent from the figure, for 
the case of periodic boundary condition, the ground state energy starts from zero at L — and grows as 
L increases. Moreover, the graph shows the existence of an inflection point around L = 5 for N = 5. In 
this region, the evaluated eigenvalues are close to the exact one {Eq — 1). Note that, the position of the 
inflection point increases for a larger set of basis functions. On the other hand, for Dirichlet boundary 
condition, the graph of the ground state energy contains a minima around L — 5 for N = 5. In this 
case, the optimal value of L (L) is the position of this minima on L axis. Also, similar to the previous 
case, L increases for larger values of N. This shows that the Rayleigh-Ritz variational formalism is only 
applicable to the case of Dirichlet boundary condition. 

Figure [3] shows the exact and optimized wave functions for N = 1,2. In fact, for N > 2 the exact and 
optimized wave functions are practically indistinguishable on the graph. The optimized ground state 
energy for N = 1 using periodic boundary condition is 1.008 (£^ xact — l) whereas this value for Dirichlet 
boundary condition are 1.136 and 1.006 for N = 1 and 2, respectively. Therefore, for this case, we need 
to utilize two oscillatory terms satisfying Dirichlet boundary condition for getting the same accuracy of 
one oscillatory term which satisfies periodic boundary condition. Moreover, the optimized wave function 



D A = E A. 



(13) 
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is also more accurate (Fig. [3]). This discrepancy is more apparent for larger values of N which will be 
addressed for the case of a double-well potential. As stated before, this is due to the fact that for the 
periodic boundary condition, the ground state wave function is not forced to vanish at x = ±L and can 
better fit itself to the exact solution which certainly is not zero there. 

Now, let us consider the situation that the potential contains an anharmonicity term (A > 0). In this 
case, the sign of k can be either positive (anharmonic oscillator) or negative (double- well potential). In 
our variational formalism, contrary to Dirichlet boundary condition (J2|), E n (L) has no minimum value 
for L > L c (Fig. [4]), where L c is the position of classical turning point. Note that, since L should be much 
larger than L c [27 [ l28 j . we consider the behavior of E n (L) only for L > L c . Figured] shows the behavior 
of the ground state energy for a double-well potential V(x) — —x 1 + x 4 and an anharmonic oscillator 
V(x) — x 2 + x 4 . As it can be seen from the figure, both plots start from zero and show a semi-flat 
region for L > L c which contains the sought after inflection point. At this point, the curvature of the 
curve vanishes. A closer look at the curve shows that this point is close to L = 3.5 for N = 8, k = 1 
and A = 1. Note that, there also exists a shallow minimum which is due to L being just large enough 
to detect the minima. Since L is too small in this case to see the overall structure of the potential, this 
minimum should be ignored. Moreover, the ground state energy of this model is Eq ~ 0.66 |26j which 
is approximately equal to the energy of the inflection point in the left part of Fig. 2J Therefore, the 
negative shallow minimum in the figure is not related to the actual value of the energy level which is 
certainly positive. 

To obtain more accurate results, we need to use more basis functions. Table Q] shows the first six 
lowest energy levels of a double- well potential (k = 1) for a small A regime (A = 0.01). We have used 
N = 100 basis functions in order to obtain the near-exact energy eigenvalues. For this case, the ground 
state energy is obtained with 82 significant digits which we have only reported it with 30 significant 
digits. As it is expected, using a fixed N, the ground state energy has the maximum precision and the 
accuracy of the higher energy levels decrease in comparison to the lower ones. To be more precise, the 
sixth eigenvalue has 79 significant digits which still is very high. Moreover, the optimum length scale 
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L for 100 basis functions is equal to 16.70762. This value strictly depends on the number of the basis 
functions and will increase for a larger set of them. Since the potential shifted by 1/(4A) is positive 
definite, some authors have reported the highly accurate positive definite eigenvalues E n + 1/(4A) |14) . 
As it can be seen from Table [T] our results are in complete agreement with them. It is worth to mention 
that, using the periodic boundary condition, the precision of 30 significant digits can be obtained using 
45 basis functions. While, for Dirichlet boundary condition, we need to use 55 basis function to get the 
same accuracy [26 . In fact, the losing of accuracy is due to the mandatory vanishing of the wave function 
at x = ±L. Note that, we have also used the optimized value of L (L) for the ground state energy to 
compute the higher energy levels E n . In fact, we have observed that L depends only on the number of 
basis functions N not on the south after energy levels. This is also true for the case of Dirichlet boundary 
condition [25] [26] |27j [28] . In principle, we need to use more basis functions to obtain the same accuracy 
for the higher energy levels. So, reduction of the accuracy of higher energy eigenvalues (using fixed N) is 
natural and is not due to the optimization procedure (Tables [T][3]). Although odd basis functions (filnix) 
have the same form for both Dirichlet and periodic boundary conditions, their optimized value of L is 
a little different which results in more accurate eigenvalues even for odd states using periodic boundary 
condition. 

For the larger values of A, the accuracy of the results increase even using the same number of basis 
functions (Table [5] and [3]), whereas the optimum value of L decreases. This is due to the decreasing 
of the classical turning points for larger A. These results are also remarkably accurate in comparison 
with the near-exact values of Refs. [HI [15]. Figure [5] shows the optimum value of L versus N for 
A = 0.01,0.03,0.1. We can also use the scaling properties of the double-well potential to find the 
properties of a general Hamiltonian H (k, A) = — p 2 — kx 2 + Ax 4 from a more simpler one which we have 
studied here H(l,(3) = —p 2 — x 2 + (3x 4 . Using the transformation of variable from x to fc 1 / 4 ^, we find 
the following scaling properties [H] [5S] 

H(k,X) = k 1/2 H(l,f3), (14) 
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Figure 2: Ground state energy of Simple Harmonic Oscillator versus L for N — 5 using periodic boundary 
condition (left) and using Dirichlet boundary condition (right). 
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Figure 3: Left, the exact and approximate ground state wave functions of SHO for L(l) = 2.52479 
and L(2) = 3.04635. Note that for N > 2 the exact and approximate wave functions are practically 
indistinguishable on the graph. Right, A^>ff =1 = \^ exact— ^ op\ versus x for periodic boundary condition 
(solid line) and Dirichlet boundary condition (dashed line). 

E(k,X) = fc 1/2 £(l,/?), (15) 
L(k,X) = fc- 1 / 4 L(l,/3), (16) 

where /? = k~ z / 2 \. Thus, it is sufficient to consider only the reduced potential (k = 1). 

4 Conclusions 

We have proposed the usage of periodic boundary condition on the basis function for accurate calculation 
of eigenvalues and eigenfunctions of a symmetric double-well potential. In this case, the Rayleigh-Ritz 
variational formalism is not applicable. Because, the graph of the energy versus the domain of the 
variable starts from zero and grows as this domain increases. We showed that this graph, because of 
the form of the potential, may contain some minima which are not related to the actual value of the 
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Figure 4: Ground state energy versus L for N = 8, A = 1, k = 1 (left), and k = — 1 (right). 
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Figure 5: L versus N for k — 1, and A = {0.01,0.03,0.1}. Circles show the position of the calculated 
inflection points and the curve is its interpolating function. 

Table 1: First six lowest energy eigenvalues of a double-well potential (k = 1) for A = 0.01, TV = 100 
and L = 16.70762. SD denotes the number of significant digits. 



n 


E n 


SD 





-23.5959513947022931175742924292 


82 


1 


-23.5959513947022931173974337194 


81 


2 


-20.8298063940006898721661249287 


80 


3 


-20.8298063940006897803867088013 


80 


4 


-18.1299111662859753878276848315 


78 


5 


-18.1299111662859531975740043181 


78 



Table 2: First six lowest energy eigenvalues of a double-well potential (k = 1) for A = 0.03, TV = 100 
and L = 13.60979. 



n 




SD 





-6.95073188927955191828148104931 


97 


1 


-6.95072754950196756189760500468 


95 


2 


-4.32728413386759375726086836212 


94 


3 


-4.32667786658379381203893295176 


94 


4 


-1.98615994840071249926930230256 


92 


5 


-1.95646376927817057309963393657 


92 
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Table 3: First six lowest energy eigenvalues of a double- well potential (k = 1) for A = 0.1, N = 100 and 
L = 11.07433. 



n E„ SD 

~Q -1.26549283721398510854595401983 lOT" 

1 -1.15305913107745006809098709688 104 

2 0.509488545436203212948452569004 103 

3 1.54354603976759862420138901373 103 

4 3.10513379668314777728015050384 101 

5 4.83611381900421025918208666909 101 

energy and should be ignored. On the other hand, we exhibit that the behavior of the energy versus the 
variable's domain shows the existence of an inflection point with vanishing curvature which is related 
to the optimum value of this domain. So, the variational procedure is finding this inflection point by 
varying the domain of the variable. We also showed that, for a fixed number of basis functions, the 
results obtained from Dirichlet boundary condition (particle-in-a-box) are less accurate in comparison 
to the usage of periodic boundary condition. In fact, by using the particle-in-a-box bases, we enforce 
the wave functions to vanish at the boundaries which results in losing the accuracy. But, the periodic 
boundary condition lets the wave functions to properly fit themselves to the exact solutions. Therefore, 
we can get the same accuracy with a smaller number of basis. For instance, for the case of k = 1, 
A = 0.01 and N = 45 basis functions, we have obtained the energy eigenvalues with 30 significant digits. 
This accuracy, for Dirichlet boundary condition, is obtainable using N — 55 particle-in-a-box basis. In 
fact, this simple and efficient method is rarely used by physicists to tackle eigenvalue problems and even 
they sometimes use some complicated and inaccurate techniques [29]. Although, we have studied only 
double-well oscillators, but the presented method is applicable quite generally for eigenvalue problems 
having polynomial potential which may correspond to physically relevant situations |30) . Moreover, this 
method can be extended and used for two- or three-dimensional cases. We expect that the priority of 
the periodic boundary condition to the Dirichlet boundary condition remains also in higher dimensions. 
In our future work, we will investigate this issue and compare our results with the exact or near-exact 
ones [31] . 
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